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Abstract 

Background: Plant root exudates have been shown to play an important role in mediating interactions between 
plant growth-promoting rhizobacteria (PGPR) and their host plants. Most investigations were performed on 
Gram-negative rhizobacteria, while much less is known about Gram-positive rhizobacteria. To elucidate early 
responses of PGPR to root exudates, we investigated changes in the transcriptome of a Gram-positive PGPR to 
plant root exudates. 

Results: Bacillus amyloliquefaciens FZB42 is a well-studied Gram-positive PGPR. To obtain a comprehensive overview 
of FZB42 gene expression in response to maize root exudates, microarray experiments were performed. A total of 
302 genes representing 8.2% of the FZB42 transcriptome showed significantly altered expression levels in the 
presence of root exudates. The majority of the genes (261) was up-regulated after incubation of FZB42 with root 
exudates, whereas only 41 genes were down-regulated. Several groups of the genes which were strongly induced 
by the root exudates are involved in metabolic pathways relating to nutrient utilization, bacterial chemotaxis and 
motility, and non-ribosomal synthesis of antimicrobial peptides and polyketides. 

Conclusions: Here we present a transcriptome analysis of the root-colonizing bacterium Bacillus amyloliquefaciens 
FZB42 in response to maize root exudates. The 302 genes identified as being differentially transcribed are proposed 
to be involved in interactions of Gram-positive bacteria with plants. 



Background 

Plant growth-promoting rhizobacteria (PGPR) are gener- 
ally referred to as a heterogeneous group of bacteria 
which colonize the rhizoplane and/or rhizosphere and 
stimulate plant growth [1,2]. PGPR have been commer- 
cially exploited as biofertilizers to improve the yield of 
crops. Some PGPR have also been successfully used as 
biocontrol agents to prevent plant diseases caused by 
phytopathogens, especially some soil-borne diseases 
[3-5]. The investigations on the interactions between 
PGPR and their host plants can not only contribute to 
our understanding of eukaryote-prokaryote relationships, 
but also have fundamental implications for designing 
new strategies to promote agricultural plant production. 
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In recent years, there is increasing evidence that plant 
root exudates play a key role in plant-microbe interac- 
tions [6-10]. Root exudates consist of an enormous 
range of compounds, among which some can attract 
beneficial associative bacteria to overcome stress situa- 
tions [8]. On the other hand, root exudates contain low 
molecular-weight carbon such as sugars and organic 
acids that primarily act as energy sources for rhizobac- 
teria and shape bacterial communities in the rhizosphere 
[11]. To date, however, it remains unclear how root exu- 
dates exert differential effects on rhizobacteria and 
which mechanisms or pathways make rhizobacteria re- 
sponsive to plant root exudates. 

Transcriptome analyses are an efficient approach to 
study host-microbe interactions at a wider scale. So far, 
the use of this approach to analyse bacterial gene expres- 
sion has been extensively used to study pathogenic 
microbes infecting their host [12]. Only a few studies 
were performed with beneficial PGPR [13-15]. Several 
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genes from Pseudomonas aeruginosa involved in metab- 
olism, chemotaxis and type II secretion were identified 
to respond to sugar-beet root exudates [13]. In another 
study, it has been suggested that the availability of par- 
ticular metabolites in root exudates, especially amino 
acids and aromatic compounds, support Pseudomonas 
putida to colonize the rhizosphere [14]. Rhizobium legu- 
minosarum was grown in the rhizospheres of its host- 
legume pea and two other non-host plants, alfalfa and 
sugar-beet. Although numerous sugar and putative com- 
plex carbohydrate transport systems are induced in the 
rhizosphere, they are less important carbon sources than 
organic acids. A common core of rhizosphere-induced 
genes was identified [15]. 

So far, studies on the impact of root exudates on 
PGPR, have been conducted with Gram-negative 
bacteria, mainly Azospirillum and Pseudomonas spp. 
[16,17]. Related studies performed with Gram-positive 
PGPR are still missing. Owing to differences in lifestyle 
and physiology, Gram-positive and Gram-negative rhizo- 
bacteria may use distinct mechanisms when interacting 
with plants. Due to their ability to produce durable 
endo-spores, bacilli are now preferred in manufacturing 
biofertilizer formulations [18], however, their successful 
application is still hampered by a lack of knowledge 
about factors determining interactions between plants 
and those bacteria, especially root colonization is far 
from being completely understood. 

Bacillus amyloliquefaciens FZB42 is a plant root- 
colonizing Gram-positive PGPR. A series of studies has 
elucidated several aspects of this rhizobacterium, particu- 
larly the molecular basis of its plant growth-promoting 
activity, which is mainly based on the production of sec- 
ondary metabolites suppressing competitive microbial 
pathogens occurring in the plant rhizosphere, the 
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secretion of the plant growth hormone auxin, and the 
synthesis of volatiles stimulating plant growth and 
induced systemic resistance (ISR) [19-21]. In the case of 
Gram-positive PGPR, however, it is still not clear how 
they maneuver their gene expression when exposed to 
plant-derived compounds. To address this question, the 
commercially established FZB42 wild type strain from 
Bacillus amyloliqufaciens was tested in this study for its 
transcriptomic responses to maize root exudates using a 
two-color DNA microarray system. 

Results and discussion 

Composition of maize root exudates 

Maize root exudates were collected from axenic hydro- 
ponic cultures and analysed by HPLC for organic acids, 
amino acids, and oligosaccharides, which have been pre- 
viously reported to be among the major ingredients in 
root exudates [8,22-24]. 

Among the compounds detected, in particular organic 
acids such as malic acid, malonic acid, succinic acid and 
trans-aconitic acid, were present at highest concentra- 
tions (Figure 1). Corroborating an earlier report [25], we 
found that lactic acid was a main constituent of maize 
root exudates. A variety of amino acids was also 
detected. Glucose and melibiose were the most promin- 
ent sugars occurring in root exudates. According to this 
analysis, most low-molecular weight organic carbon 
appeared to be present in the form of organic acids. 

Overall changes in gene expression in response to root 
exudates 

In the rhizosphere, root exudates may occur at high con- 
centrations in certain microenvironments, e.g. in vicinity 
of root tips [26], but their concentration in specific 
niches of the environment is unknown. Therefore, the 
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Figure 1 Composition and concentration of the maize root exudates. Exudates collected from the roots of maize seedlings were assayed by 
HPLC. Organic acids, amino acids, and carbohydrates were quantified and those which had a concentration of >0.1 umol g-1 (dry weight) were 
included in the graph. Proline, a known constituent of maize root exudate, was not detected since the derivatization reagent (OPA) used reacts 
only with primary amino groups. 
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choice of a physiologically relevant concentration of exu- 
dates to be used for microarray experiments can only be 
tentative. Based on a previous study on changes in the 
proteomics of FZB42 [27], three exudate concentrations 
(0.25 gl" 1 , 0.5 gl" 1 and 1.0 gl" 1 ) were applied to liquid 
cultures of FZB42, and bacterial cells were harvested for 
RNA extraction at two growth stages (OD 600 = 1.0 and 
OD 60 o = 3.0). For simplicity, the two population densities 
were referred to as OD1.0 and OD3.0 throughout this 
paper, respectively. A concentration of 0.25 g l" 1 was suf- 
ficient to result in a significant response of FZB42 tran- 
scriptome. When bacteria were cultured at OD3.0 the 
number of up-regulated genes gradually decreased with 
increasing root exudate concentration, suggesting that 
some compounds need to occur at lower abundance to 
induce gene expression, or that gene transcription in 
general may be suppressed at high concentrations of 
some exudates components (Figure 2). More transcripts 
were significantly altered (q<0.01) at the transition to 
the stationary growth (OD 3.0) than at the exponential 
growth (OD1.0) (Figure 2), suggesting that OD 3.0 was a 
sampling point which reflected more clearly the effect of 
root exudates on FZB42 than OD1.0. For these reasons, 
the exudate concentration of 0.25 gl" 1 and the OD3.0 
for harvesting of cells were used for all subsequent 
microarray experiments. 

Six independent experiments were performed and the 
genes whose transcription fulfilled the condition of 
yielding a q value not greater than 0.01 (q < 0.01) and a 



fold change not less than 1.5 (FCH > 1.5) were regarded 
as being significantly influenced by root exudates. A 
total of 302 genes, representing 8.2% of the FZB42 tran- 
scriptome, were significantly regulated in their transcript 
levels by the applied root exudates (see supplemental 
material Additional file 1: Tables SI, S2, and S3). The 
majority of these genes (261 genes) was up-regulated, 
whereas only 41 genes were down- regulated (Figure 3). 
Although most of the regulated genes have been func- 
tionally annotated, a significant proportion (-23%) 
remained of unknown function, among which 19 genes 
were unique for FZB42. In addition, 44 genes (-15%) 
encoded either hypothetical proteins or proteins with 
putative functions (Figure 3). The distribution in various 
functional categories of all the gene with known (189 
genes) or putative (44 genes) products are summarized 
in Figure 4. 



Validation of microarray result by real-time PCR 

Nine up-regulated genes with different levels of fold 
changes in expression (1.5 - 5.2 fold) were chosen to be 
evaluated by quantitative real-time PCR. All these genes 
were confirmed to be significantly up-regulated in the 
presence of root exudates (Figure 5). The fold change of 
each gene revealed by real-time PCR was similar to that 
obtained in the microarray experiments (Figure 5). In 
summary, the real-time PCR suggested that the micro- 
array data were reliable. 




0 >1. 0-fold up-regulated 
D >1. 0-fold down-regulated 



0 >1. 5-fold up-regulated 
D >1. 5-fold down-regulated 



Figure 2 Number of FZB42 genes altered in transcription in response to root exudates at different exudate concentrations and cell 
densities. Maize root exudates were supplemented in three concentrations (0.25 mg/ml, 0.5 mg/ml and 1.0 mg/ml) to FZB42 cultures and total 
RNA was prepared from the bacterial cells harvested at two optical densities (OD 600 = 1.0 and OD 600 = 3.0). Genes significantly altered in 
transcription (q < 0.01 and fold change >1.5) by presence of root exudates are represented in the figure. 
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Figure 3 Overview of groups of the 302 genes altered in 
transcription by root exudates. A total of 302 genes were 
significantly altered (q < 0.01 and fold change >1.5) in transcription 
by the maize root exudates. "Up" indicates genes that were 
up-regulated in presence root exudates, while "down" the ones that 
were down-regulated by the root exudates. The genes encoding a 
product with known or unknown function and those encoding a 
hypothetical protein were indicated. The number of genes of each 
section and their percentage is depicted. 



The regulated genes with known function 

Among the 302 genes with significantly altered expres- 
sion by root exudates, 189 were annotated with known 
functions. These were categorized into various classes 
[28], such as cell envelope and cellular processes, inter- 
mediary metabolism, information pathway and other 
functions . In these categories, three groups (Table 1) 
contained the largest numbers of genes and at least one 
third of the genes within these groups had a fold change 
of >2.0. This suggests that these three groups of genes 
were strongly affected by root exudates: 

i) The transcription of 46 genes involved in carbon 
and nitrogen utilization was altered in response to root 
exudates, with 43 of them being up-regulated. These 46 
genes were involved in different aspects of the metabol- 
ism of carbohydrates, amino acids and related metabo- 
lites. To obtain a more comprehensive understanding of 
their relevance in the metabolic context, the genes were 
mapped into the KEGG pathway and a representation of 
metabolic pathways was constructed (Figure 6). A total 
of 12 genes encoding enzymes involved in the Embden- 
Meyerhof-Parnas (EMP) pathway (including pgi encod- 
ing for glucose-6-phosphate isomerase) and the TCA 
cycle were significantly up-regulated. These genes cov- 
ered almost the entire glycolysis and TCA pathway. 



Nearly a quarter of the genes with altered transcription 
(46 out of 189) were involved in uptake or utilization of 
nutrients. This observation corroborated that root exu- 
dates serve as energy sources in the interaction between 
roots and rhizobacteria. 

Among the up-regulated genes, glvA, glvC and glvR 
showed the highest fold change (glvA: 5.2-fold up-regu- 
lated, glvC: 2.5-fold up-regulated, glvR: 4.4-fold up-regu- 
lated). The enhancement of glvA expression was also 
validated by real-time PCR as well as by the proteomic 
data (unpublished). The three genes comprise the glv op- 
eron (glvA-glvR-glvC), which is responsible for maltose 
dissimilation and positively regulated by maltose [29]. 
The significant up-regulation of these genes indicated 
that maltose was present in the exudates, which was 
confirmed by the HPLC analysis (Figure 1). 

The genes involved in inositol metabolism (iolA, iolB, 
iolC, iolD, iolE, iolF, iolG, ioll, MS) were also up-regu- 
lated, mainly with a fold change of >2.0 (Figure 6). Ex- 
cept iolS, which is involved in the regulation of inositol 
catabolism, the other eight genes are members of the iol 
operon. The increased transcription of iolA and iolD was 
further confirmed by real-time PCR whereas the en- 
hancement of iolB and ML was validated by a proteo- 
mics approach (unpublished data). The activation of 
nine genes indicated the presence of inositol in the exu- 
dates, which has also been verified by HPLC. 

ii) A second group of genes with a higher fold change 
were those associated with sensing, chemotaxis, motility 
and biofilm formation (Table 2). These processes are 
crucial for bacterial colonization of roots. The recogni- 
tion of signals released from roots and rhizobacteria is 
the first step of the establishment of a mutual cross -talk 
[30]. Once plant signals have been perceived, bacteria 
move towards the plant root to establish in the rhizo- 
sphere [31-34]. Bacterial motility in the rhizosphere 
involves several processes such as chemotaxis, flagella- 
driven motility, swarming, and production of surfactants 
[35-38]. The observed transcriptional changes of genes 
required for chemotaxis {cheC, cheD) and motility {hag 
fliD, fliP and flgM) indicated that root exudates contain 
compounds that induce attraction of FZB42 cells to 
roots. 

Biofilm formation has been documented to be involved 
in directing or modulating efficient colonization by 
PGPR [39,40]. Biofilms can also provide the plant root 
system with a protective barrier against attack of patho- 
genic microbes [35]. Two B. amyloliquefaciens genes 
involved in biofilm formation, ycmA and luxS, were 
enhanced by maize root exudates (Table 2, Additional 
file 1: Table SI). The gene luxS, required for synthesis 
of the quorum-sensing signaling molecule autoinducer- 
2 (AI-2) [41], is involved in biofilm formation of 
pathogenic Streptococcus species [42-44] and the 
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Figure 4 Distribution in various functional categories of the genes altered in transcription by root exudates. Among the 302 genes 
altered in transcription by maize root exudates at OD3.0, those with known (189 genes) or putative (44 genes) products were classified according 
to their function. The percentage of each group is indicated. 

V J 



probiotic B. subtilis natto [45]. The gene ycmA has 
also been indentified to be involved in facilitating bio- 
film formation [46,47]. The increased transcription of 
luxS and ycmA indicated that biofilm formation of 
FZB42 could be enhanced by some compounds 
present in root exudates. 



iii) The third functional group with the highest num- 
ber of genes induced by root exudates was associated 
with the non-ribosomal synthesis of secondary metabo- 
lites with antimicrobial action (Table 3). Producing sec- 
ondary metabolites suppressing deleterious microbes in 
the rhizosphere is an established mechanism of 
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Figure 5 Fold-change of differentially expressed genes selected for validation by Real-time PCR. The fold changes revealed by real-time 
PCR of the selected genes were determined using the software REST. Three repeats were performed for each gene. For comparison, the fold 
changes obtained in microarray analysis were shown in parenthesis below each specific gene. The boxes represent the distance between the 
25th and the 75th percentile. The lines in the boxes represent the median gene expression. Whiskers represent the minimum and maximum 
observations. 
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Table 1 Functional categories* of the FZB42 genes 
significantly regulated by the maize root exudates and 



with known functions 

Classification code_Functional category Nr. of the 

genes included 

1 cell envelope and cellular processes 58 

1.7_ Cell division 6 

1.1_ Cell wall 5 

1.4_ Membrane bioenergetics 7 

1.5_ Mobility and chemotaxis 6 

1.3_ Sensors (signal transduction) 2 

1.6_ Protein secretion 5 

1.8_ Sporulation 7 

1 .1_ Transformation/competence 2 

1.2_ Transport/binding proteins and lipoproteins 18 

2_intermediary metabolism 59 

2.1_Metabolism of carbohydrates and 34 
related molecules 

2.2_ Metabolism of amino acids and 12 
related molecules 

2.5_ Metabolism of coenzymes and prosthetic groups 4 

2.4_ Metabolism of lipids 5 

2.3_ Metabolism of nucleotides and nucleic acids 4 

3_information pathways 45 

3.3_ DNA recombination 1 

3.1_ DNA replication 3 

3.8_ Protein modification 2 

3.7_ Protein synthesis 20 

3.6_ RNA modification 1 

3.5_ RNA synthesis 18 

4_other functions 27 

4.1 _ Adaptation to atypical conditions 6 

4.2_ Detoxification 4 

4.6_ Miscellaneous 3 

4.4_ Phage-related functions 1 

4.3_ Antibiotic production 13 

In total 189 



The categories in which more than one third of the genes had a fold change 
of >2.0 were highlighted in bold text (Refer to experiment "Response to RE": 
E-MEXP-3421). *The genes were categorized according to [28]. 



biocontrol adopted by B. amyloliquefaciens FZB42 on 
plants [19,48,49]. The majority of the induced genes are 
devoted to the synthesis of two polyketide antibiotics, 
bacillaene and difficidin. Some components in the exu- 
dates could stimulate the production of these two anti- 
biotics, which have been demonstrated to be able to 
protect orchard trees from fire blight disease caused by 
Erwinia amylovora [49]. 

Another two genes, mlnH and fenE, were also in- 
duced, which are known to participate in non-ribosmal 



biosynthesis of macrolactin and fengycin, respectively. 
Macrolactin, a polyketide product found in FZB42, 
has activity against some Gram-positive bacteria [50], 
while fengycin can act against phytopathogenic fungi 
in a synergistic manner with bacillomycin D [19,51]. 

In addition, two genes encoding surfactin synthe- 
tase were also activated by root exudates (Table 3). 
Surfactin is one of Bacillus cyclic lipopeptides, dis- 
playing antiviral and antibacterial activities. In Arabi- 
dopsis it has been shown that the ability of Bacillus 
to synthesize surfactin can reduce the invasion of 
Pseudomonas syringae [30]. although it is not yet 
clear whether the protective effect resulted directly 
from the antibacterial activity of surfactin or from its 
biofilm-related properties. Surfactin is crucially 
involved in the motility of Bacillus by reducing sur- 
face tensions [36,37,52] and contributing to biofilm 
formation on Arabidopsis roots [30]. It has also been 
demonstrated that surfactin production of FZB42 was 
enhanced when colonizing the duckweed plant Lemna 
minor [21]. It can be expected that up-regulation of 
srfAC and srfAD may contribute to the protective 
role of surfactin against plant pathogens. 

The regulated genes with putative function 

Among the 302 genes significantly altered in tran- 
scription by root exudates, 44 were annotated to en- 
code a putative enzyme or a hypothetical protein. 
Similar to the genes with known function, these 44 
genes fell into three categories: metabolism of carbo- 
hydrates and related molecules, metabolism of amino 
acids and related molecules, and transport/binding 
proteins and lipoproteins (Additional file 1: Table S2). 
Some of the 44 genes were closely associated with 
plant-microbe interactions. For example, the tran- 
scription of ydjL, nowadays renamed bdhA, encoding 
acetoin reductase/butanediol dehydrogenase [53], was 
1.5-fold enhanced by root exudates. 2, 3-Butanediol is 
a volatile organic compound released by PGPR and 
able to promote significantly plant growth [54]. The 
expression of the gene epsE, residing in a 15-gene 
operon epsA-O, was also enhanced by root exudates. 
EpsE is involved in formation of biofilm by arresting 
flagellar rotation of cells embedded in biofilm matrix 
[55]. Another activated gene was dfhY, which encodes a 
hypothetical protein. Like other induced genes known to 
be involved in antibiotic production such as dfhF, dfnG, 
dfnl and dfnj (Table 3), dfhY is part of the gene cluster re- 
sponsible for synthesis of the polyketide antibiotic diffici- 
din. It is worth mentioning that antibiotic production is 
energetically very costly and its strict control is a clear evo- 
lutionary advantage. 

In contrast to a few genes significantly altered during 
the exponential phase (OD1.0), hundreds of genes were 
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Figure 6 A subset of the up-regulated genes with known function in response to maize root exudates. The significantly up-regulated 
genes by the root exudates were mapped in the KEGG pathway and the diagram was accordingly adapted. The products encoded by the 
up-regulated genes were highlighted in red, whilst the down-regulated YadH was highlighted in green. CM stands for cell membrane. 



differentially expressed in presence of root exudates dur- 
ing transition to stationary growth phase (OD3.0). Such 
a difference may not be surprising. The transcription of 
most bacterial genes during the exponential growth 
phase is typically initiated by RNA polymerase holoen- 
zyme carrying the housekeeping transcription factor cA 
while in the stationary phase, transcription is mainly 
accomplished by RNAP carrying alternative sigma 



Table 2 FZB42 genes significantly induced by maize root 
exudates and involved in mobility and chemotaxis (Refer 
to experiment "Response to RE": E-MEXP-3421) 



Gene 


Fold change 


Classification code_function involved 


MM 


2.0 


1.5_ Mobility and chemotaxis 


fliP 


1.7 


1.5_ Mobility and chemotaxis 


cheC 


1.7 


1.5_ Mobility and chemotaxis 


cheD 


-1.5 


1.5_ Mobility and chemotaxis 


hag 


3.6 


1.5_ Mobility and chemotaxis 


flgM 


1.7 


1.5_ Mobility and chemotaxis 


luxS 


1.7 


1.3_ Sensors (signal transduction) 


ymcA 


2.5 


1.3_ Sensors (signal transduction) 



Table 3 FZB42 genes which were significantly induced by 
maize root exudates and involved in antibiotic 
production (Refer to experiment "Response to RE": 
E-MEXP-3421) 



Gene 


Product 


Fold 
change 


boeE 


malonyl-CoA-[acyl-carrier protein] transacylase BaeE 


1.6 


boel 


enoyl-CoA-hydratase Bael 


2.2 


baeL 


polyketide synthase BaeL 


1.9 


baeN 


hybrid NRPS/PKS BaeN 


1.5 


baeR 


polyketide synthase BaeR 


2.3 


dfnJ 


modular polyketide synthase of type 1 DfnJ 


2 


dfnl 


modular polyketide synthase of type 1 Dfnl 


1.7 


dfnG 


modular polyketide synthase of type 1 DfnG 


2 


dfnF 


modular polyketide synthase of type 1 DfnF 


2.4 


mlnH 


polyketide synthase of type 1 MlnH 


1.5 


fenE 


fengycin synthetase FenE 


1.5 


srfAD 


surfactin synthetase D SrfAD 


1.9 



srfAC surfactin synthetase C SrfAC 1.7 
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factors allowing to adapt to a permanently changing en- 
vironment. The extracytoplasmic-function (ECF) sigma 
factor W was enhanced in presence of root-exudate 
(Figure 5). SigW is known as being expressed in early 
stationary growth-phase and induced by various cell wall 
antibiotics, alkaline shock, and other stresses affecting 
the cell envelope. It controls a large "antibiosis" regulon 
involved in mediating resistance to various antibiotics 
including fosfomycin and the antibiotic peptides sublan- 
cin and SdpC [56]. It has been observed that many 
virulence-associated factors influence the colonization, 
persistence and spreading mechanisms of the human 
pathogen Streptococcus pyogenes in a growth phase- 
dependent manner [57-59]. Likewise, rhizobacteria may 
employ an early stationary phase-related mechanism to 
favor expression of those genes that mediate rhizosphere 
competence. 

Effect of soil extract 

To simulate in part the conditions that bacteria experi- 
ence in the rhizosphere soil, 10% soil extract was added 
to the culture media. Additional microarray experiments 
were performed in a similar way as before to investigate 
the effect of the soil extract on gene expression of 
FZB42. The result showed that no gene was significantly 
up-regulated by the soil extract during exponential 
growth phase of OD1.0, whereas five genes were 
repressed in the presence of the soil extract at OD3.0 
(Table 4). This negligible number of genes that were dif- 
ferentially transcribed indicates that the supplement of a 
soil extract did not have major effects on gene transcrip- 
tion under the growth conditions used. 

Effect of exudates prepared from maize plants colonized 
by FZB42 

Typically, most root exudates studied were collected 
from plants grown in axenic systems. The release of root 
exudates is not only determined by the plant species, but 
also by plant age, physiological status, and the biotic en- 
vironment that plants thrive including the rhizosphere 
microflora that influence the composition and quantity of 
root exudates [60-66]. It was reported that P. aeruginosa 
produces A/-acyl homoserine lactone (AHL) signaling 



compounds that induce changes in the root exudation of 
Medicago truncatula [67]. Exudate compounds that are 
specifically induced or repressed by rhizobacteria may in 
turn affect bacterial gene expression. Such an effect can- 
not be demonstrated using root exudates collected from 
a gnotobiotic system, therefore, a batch of "interaction 
exudates (IE)" was collected from maize roots which were 
previously inoculated with FZB42. 

The transcriptional responses of FZB42 to the IE were 
compared with responses to the root exudates (RE) col- 
lected from axenic culture. No significant differences 
(q<0.01 and FCH>1.5) were found between the effect 
of IE and RE at OD1.0, while four genes were differen- 
tially expressed at OD3.0 (Additional file 2: Table S5). 
When a less stringent selection filter was applied 
(q<0.05 and FCH> 1.5), a total of nine genes were dif- 
ferentially expressed (Additional file 2: Table S5). The 
four genes, significantly enhanced in presence of FZB42 
at maize roots, encode enzymes involved in the degrad- 
ation of macromolecules or cellular compounds, such as 
ggt, nprE, clpP, RBAM00438 (ycsN). Among all four 
genes, expression of the ggt gene was found most 
enhanced, bearing a fold change of 2.2 in presence of 
the rhizobacterium (Additional file 2: Table S5). GGT, y- 
glutamyltranspeptidase (GGT) (EC 2.3.2.2) catalyzes the 
hydrolysis of y-glutamyl compounds, such as glutathione 
(GSH), and the transfer of y-glutamyl moieties to amino 
acids and peptides. The nprE gene, which is mainly 
expressed during early stationary phase, encodes extra- 
cellular neutral protease involved in degradation of pro- 
teins and peptides. The peptidase ClpP, encoded by the 
clpP gene, can associate with the ATPases ClpC, ClpE, 
and ClpX, thereby forming a substrate specific channel 
for several regulatory proteins directing spore formation 
or genetic competence in bacilli. RBAM00438 is a mem- 
ber of the aldo-keto reductases (AKRs) superfamily of 
soluble NAD(P)(H) oxidoreductases whose chief purpose 
is to reduce aldehydes and ketones to primary and sec- 
ondary alcohols. At present, it remains questionable if 
those gene products are linked with any specific process 
triggered by the IE. The number of the genes obtained 
was much less than expected. We conclude that possible 
differences between the transcriptome responses to these 



Table 4 FZB42 genes repressed by soil extract at OD3.0 (Refer to experiment "Response to SE": E-MEXP-3551) 



Gene 


Fold change 


Product 


Function involved 


ypeO 


-2.6 


hypothetical protein YpeQ 


unknown 


yurV 


-2.4 


iron-sulfur cofactor synthesis protein nifU homolog YurV 


miscellaneous 


iolS 


-2.2 


inositol utilization protein S (IolS) 


metabolism of carbohydrates and related molecules 


yooA 


-2.0 


conserved hypothetical protein YaaA 


unknown 


ahpF 


-2.0 


alkyl hydroperoxide reductase (large subunit) and NADH 
dehydrogenase AhpF 


detoxification 
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two exudate samples are either very rare or too subtle to 
be revealed sufficiently by two-color microarrays. 

One drawback of the present investigation is that some 
effects of the root exudates may have been masked by 
components of the 1 C medium and therefore did not 
result in altered gene expression. On the other hand, 
using 0.25 mg exudates per ml medium, some compo- 
nents in the exudates may have been diluted to a level at 
which they no longer show detectable effect on bacterial 
gene expression. It has been reported that the rhizo- 
sphere is a very heterogeneous soil volume, with some 
regions being "hotspots" of root exudation and bacterial 
colonization. In natural environments, bacterial popula- 
tions are likely to be exposed to different concentration 
of exudates along the root axis [68,69]. 

It needs to be mentioned that the exudates used in this 
study were a pooled mixture of the samples collected 
within seven days from maize roots (see Methods). It 
has not yet been described to which extent the compos- 
ition of root exudates is affected by the developmental 
stage of a plant [70] and therefore the presented bacter- 
ial responses cannot be assigned to a particular physio- 
logical state of the host plant. This question may be 
addressed by performing bacterial transcriptome ana- 
lyses in response to exudates collected at different time 
points during plant development. Such an approach may 
be helpful to elucidate the progression of the plant- 
bacteria association during the plant development. 

In summary, this microarray work reflects the interac- 
tions between a Gram-positive rhizobacterium and its 
host plant in a genome-scale perspective. Critical target 
genes and pathways for further investigations of the 
interaction were identified. Given the limited reports on 
transcriptomic analysis of rhizobacteria in response to 
their host plants [13-15], the results provided a valuable 
insight into PGPR behaviour in the rhizosphere. About 
10% of the total number of genes were found up- 
regulated in presence of root exudate during transition 
to stationary growth phase. In addition to the findings 
corroborating previous transcriptome analyses per- 
formed in Gram-negative bacteria, we could demon- 
strate that presence of root exudate induced expression 
of numerous genes involved in non-ribosomal synthesis 
of secondary metabolites with antifungal and antibacter- 
ial action. We hypothesize that competitive colonization 
at plant root surfaces by FZB42 might be supported by 
enhanced synthesis of antimicrobial compounds. 

Conclusions 

Using the data from six independent micro array experi- 
ments, differentially transcribed genes of the PGPR B. 
amyloliquefaciens FZB42 were identified and their 
known or putative functions were related to their asso- 
ciative behavior with regard to interactions with maize 



roots. A large group of genes specifically expressed sug- 
gested that root exudates serve primarily as a source of 
carbon and energy for FZB42. Another group of genes 
significantly induced by plant root exudates encode the 
non-ribosomal synthesis of antimicrobial secondary 
metabolites. It is possible that enhanced synthesis of 
antimicrobial compounds might suppress the competing 
phytopathogenic organisms growing within the plant 
rhizosphere. However, direct evidence for occurrence of 
those compounds in vicinity of plant rhizosphere 
remains to be accomplished. The addition of soil 
extracts to the growth medium showed no major effect 
on gene expression of FZB42. Similarly, the results 
obtained with the "interaction exudates" collected from 
the maize roots inoculated with FZB42 did not indicate 
altered effects on gene expression compared with that of 
common root exudates collected in the gnotobiotic 
system. 

Methods 

Root exudates collection and analysis 

Maize seeds (Saaten-Union, Germany) were surface- 
sterilized and germinated as described previously [21]. 
Root exudates were collected from the maize seedlings 
grown in an axenic system with sterile water (1:1 dis- 
tilled water and tap water, v/v). Forty germinated seeds 
harboring a main root of at least 2 cm length were trans- 
ferred into test tubes filled with 2 ml of autoclaved 
water, with the maize seeds being placed just above the 
water surface. The tubes were kept under sterile condi- 
tions and maintained in a plant growth room (16-h 
light/8-h dark) at 24°C for 8 days. In the first two days, 
water was supplemented to the tubes, and seedlings 
were pulled to a higher position to ensure that the maize 
seeds were always above the water surface as the roots 
elongated. From the third day on, the water containing 
the exudates was collected and the tubes were refilled 
with sterile water. Sampling was performed every day 
until the eighth day after transferring the seedlings. Each 
collection were kept separate, from which a 100 \iL ali- 
quot was taken and spread on a solid LB media to check 
for contamination. The contaminated samples were 
discarded. 

To collect the "interaction exudates (IE)", the germi- 
nated maize seeds were inoculated with FZB42 as 
described previously [21] before transferring the test 
tubes. Afterward the maize was grown and the exudates 
were prepared in the same way as described above. 

The collected exudates were pooled, freeze-dried and 
stored at -20°C. Before use, the lyophilized exudates 
were weighted, and dissolved in a certain volume of dis- 
tilled water. The obtained exudates solution was centri- 
fuged to remove any insoluble constituents. The 
supernatant was filter-sterilized and the resulting stock 
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exudates were stored in dark at -80°C The final concen- 
tration of the exudates in the culture vessel was gener- 
ally adjusted to 0.25 gL 1 . Chemical analysis of the root 
exudates was performed as described previously [71]: 
amino acids were determined using a Shimadzu HPLC 
system. 40 \iL samples were derivatized with 160 (il OPA 
(o-phthaldialdehyde) reagent and 20 \iL of the resulting 
mixture were injected and separated on a GROM-SIL 
OPA-3 column using solvent gradient elution by solvent 
A (25 mM phosphate buffer pH 7.2 with 0.75% tetra- 
hydrofuran) and solvent B (methanol to acetonitrile to 
25 mM phosphate buffer 7.2 [35 : 15 : 50/v : v : v]). Gra- 
dient profile: 0-2 min, 0% B; 2-10 min, 0%-50% B; 10- 
15 min, 50-60% B; 15-20 min, 60-100% B; 20-25 min, 
100% B; 25-26 min, 100%-0% B; 26-35 min, 0% B. The 
flow-rate was 1 mL min" 1 . Subsequent fluorescence de- 
tection of the derivatives was performed at an excitation 
wavelength of 330 nm and 450 nm. Organic acids were 
determined by means of ion chromatography (Dionex 
IonPac AS 11 HC column) using a gradient ranging 
from 4 mM to 80 mM KOH. Organic acids were identi- 
fied by comparison of retention time with known stan- 
dards. Sugars were determined by GC-TOF-MS. A 
lyophilized 75 \xL aliquot of root exudates was dissolved 
in 50 mL methoxyamine hydrochloride in dry pyrididine 
and derivatized for 2 h at 37°C followed by 30 min. treat- 
ment with 50 \iL N-methyl-N-trifluoroacetamide at 37° 
C. A volume of 1 \iL was injected into the GC column. 

Microarray design 

The Bam4kOLI microarray was designed based on the 
sequenced complete genome of B. amyloliquefaciens 
FZB42 [27] (Additional file 3: Table S6). The array con- 
tained 3931 50-70mer oligonucleotides representing 
predicted protein-encoding genes and a set of small 
non-coding RNA genes of FZB42. In addition, the array 
included stringency controls with 71%, 80% and 89% 
identity to the native sequences of five genes, dnaA, 
rpsL, rpsO, rpsP, and rprnl, to monitor the extent of 
cross hybridization. The array also contained alien DNA 
oligonucleotides for four antibiotic resistance genes 
(Ern r , Cm r , Nm r and Spc r ) and eight spiking controls as 
well as one empty control (nothing spotted). All oligo- 
nucleotide probes were printed in four replicates. Micro- 
arrays were produced and processed as described 
previously [72]. 

Oligonucleotides were designed using the Oligo De- 
signer software (Bioinformatics Resource Facility, CeBi- 
Tec, Bielefeld University). Melting temperatures of the 
oligonucleotides were calculated based on %GC and 
oligo length, ranging from 73°C to 83°C (optimal 78°C). 
Salt concentration was set to 0.1 M. QGramMatch was 
used to analyse uniqueness of the oligos. 



Experimental design 

The experiment designs of FZB42 in response to various 
conditions are summarized in Additional file 3: Table S6. 
Independent experiments were used as biological repli- 
cates. In all comparisons dye-swap were carried out to 
minimize the effect of dye biases. 

1 C medium (0.7% w/v pancreatic digest of casein, 
0.3% w/v papain digest of soya flour, 0.5% w/v NaCl) 
containing 0.1% glucose was used in all experiments. Ex- 
cept the controls of the experiment "Response to SE" 
(Additional file 3: Table S6), 10% soil extract was also 
supplemented in the media. Soil extract was prepared by 
extracting 500 g dried, fertile garden soil with one litre 
distilled water for 2 hrs and autoclaving. After cooling 
down, the supernatant was filtered with 0.22 \im Nucle- 
pore unit and then stored at 4°C until use. 

Total RNA preparation 

One overnight colony of FZB42 was inoculated into 1 C 
medium plus 0.1% glucose and then shaken at 210 rpm 
at 24°C. After 14 hours the obtained preculture was used 
to inoculate a new 1 C medium (containing 0.1% glu- 
cose) plus the corresponding solution to be studied 
(maize root exudates, soil extract, or interaction exu- 
dates. See Additional file 3: Table S6). The main cultures 
were grown at 24°C until they reached late exponential 
growth phase (OD 1.0) and/or the transition to station- 
ary phase (OD3.0, see Additional file 4: Figure SI). 

The FZB42 cells of OD1.0 or OD3.0 were harvested 
for preparation of total RNA. A volume of 15 ml of the 
culture was mixed with 7.5 ml "killing buffer" (20 mM 
Tris-HCl, 5 mM MgCl2, 20 mM NaN3, pH 7.5) and 
then centrifuged at 5,000 rpm for 3 minutes at room 
temperature. The pellet was washed once more with 
1 ml "killing buffer" and then immediately frozen in li- 
quid nitrogen. The frozen cell pellets were stored at 
-80°C until RNA isolation. 

Isolation of RNA was performed using the Nucleo 
Spin® RNA L (Macherey Nagel) according to the manu- 
facturers instructions. The isolated RNA was addition- 
ally digested with DNasel to avoid possible trace DNA 
contamination. After ethanol precipitation RNA pellets 
were resuspended in 300 \A RNase-free water. The con- 
centration of total RNA was spectrophotometrically 
determined, whereas its quality was checked on a 1.5% 
RNA agarose gel in 1 x MEN buffer (20 mM MOPS; 
1 mM EDTA, 5 mM NaAc; pH7.0) with 16% 
formaldehyde. 

Synthesis of labeled cDNA, hybridization and image 
acquisition 

Synthesis of first-strand cDNA, microarray hybridization 
and image acquisition were performed in CeBiTec, the 
Center for Biotechnology at Bielefeld University. Briefly, 
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aminoallyl-modified first-strand cDNAs were synthe- 
sized by reverse transcription according to DeRisi et al 
[73]. and then coupled with Cy3- and Cy5-A/-hydroxy- 
succinimidyl ester dyes (GE Healthcare, Little Chalfont, 
UK). After hybridization using the HS4800 hybridization 
station (Tecan Trading AG, Switzerland), slides were 
scanned with a pixel size of 10 \im using the LS 
Reloaded microarray scanner (Tecan Trading AG, 
Switzerland). 

Data processing 

The microarray data obtained was analysed by using the 
EMMA 2.8.2 software [74]. The mean signal intensity 
(Ai) was calculated for each spot using the formula A- x = 
\og 2 (RiGi) 0 - 5 [75]. #i = 4hi(i)-£&hi(i) and Gi = / ch2(i) - 
Bgcb2(i)> where 7 c hi(i) or 7 c h2(i) is the intensity of a spot in 
channel 1 or channel 2, and £gchi(i) or Bg ch2 (i) is the 
background intensity of a spot in channell or channel 2, 
respectively. The log 2 value of the ratio of signal inten- 
sities (Mi) was calculated for each spot using the for- 
mula Mi = log 2 C#i/Gi). Spots were flagged as "empty" if 
R < 0.5 in both channels, where R = (signal mean-back- 
ground mean) /background standard deviation [76]. The 
raw data were normalized by the method of LOWESS 
(locally weighted scattered plot smoothing). A signifi- 
cance test was performed by the method of false discov- 
ery rate (FDR) control and the adjusted p-value defined 
by FDR was called q- value [77,78]. 

An arbitrary cutoff, fold change (FCH) greater than 
1.5, was applied to the genes with a q-value of <0.01. 
Only those genes which meet both filter conditions 
(q<0.01 & FCH > 1.5) were regarded to be significantly 
differentially expressed. 

Real-time PCR 

The first-strand cDNA was obtained by reverse tran- 
scription with RevertAid™ Premium Reverse Tran- 
scriptase (Fermentas, St. Leon-Rot, Germany), using 
random hexamers as primers. Oligonucleotide primers 
were designed by the software PrimerExpress and listed 
in supplemental materials (Additional files 1: Table S4). 
Real-time PCR was performed with SYBR® Green PCR 
Master Mix kit (Carlsbad, California, USA) using 7500 
Fast Real-Time PCR System (Carlsbad, California, USA) 
according to the manufacturers' instructions. As an in- 
ternal control, the housekeeping gene gyrA was used as 
its expression was not significantly altered in all micro- 
array experiments. Three technical replicates were car- 
ried out for each target gene. Quantification was 
analysed based on the threshold cycle (Ct) values as 
described by Pfaffl [79]. 

The raw data of the Micro-array experiments, 
described here, are available in the ArrayExpress data- 
base under the accession numbers: E-MEXP-3421, E- 



MEXP-3550, E-MEXP-3551, E-MEXP-3553, E-MEXP- 
3554, respectively (see also Additional file 3: Table S6). 
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